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. We present a systematic exposition of the Lagrangian field theory for the mas- 

"q \ sive spin-two field generated in higher-derivative gravity upon reduction to a 

O^' second-order theory by means of the appropriate Legendre transformation. It 

. has been noticed by various authors that this nonlinear field overcomes the well 

known inconsistency of the theory for a linear massive spin-two field interacting 
with Einstein's gravity. Starting from a Lagrangian quadratically depending on 
^ . the Ricci tensor of the metric, we explore the two possible second-order pictures 

usually called "(Helmholtz-) Jordan frame" and "Einstein frame". In spite of their 
mathematical equivalence, the two frames have different structural properties: in 
Einstein frame, the spin-two field is minimally coupled to gravity, while in the 
other frame it is necessarily coupled to the curvature, without a separate kinetic 
term. We prove that the theory admits a unique and linearly stable ground state 
solution, and that the equations of motion are consistent, showing that these 
results can be obtained independently in either frame (each frame therefore pro- 
vides a self-contained theory). The full equations of motion and the (variational) 
energy-momentum tensor for the spin-two field in Einstein frame are given, and a 
simple but nontrivial exact solution to these equations is found. The comparison 
of the energy-momentum tensors for the spin-two field in the two frames suggests 
that the Einstein frame is physically more acceptable. We point out that the 
energy-momentum tensor generated by the Lagrangian of the linearized theory is 
unrelated to the corresponding tensor of the full theory. It is then argued that 
the ghost-like nature of the nonlinear spin-two field, found long ago in the linear 
approximation, may not be so harmful to classical stability issues, as has been 
expected. 
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1 Introduction 



A consistent theory of a gravitationally interacting spin-two field could not be developed 
until a significant progress was made in an apparently unrelated subject, i.e. higher- 
derivative metric theories of gravity. It is well known that a single linear spin-two field 
cannot be consistently coupled to gravity. It is therefore a common belief that Nature 
avoids the consistency problem by simply not creating fundamental spin-two (nor higher 
spin) fields except gravity itself. Nevertheless the subject has remained fascinating over 
decades and some authors have studied various aspects of linear spin-two fields ^ |^ , 
in particular their dynamics in Einstein spaces 0, |^. 

On the other hand, higher-derivative metric theories of gravity, where the Lagrangian 
is a scalar nonlinear function of the curvature tensor (hence in this paper they are 
named nonlinear gravity theories, NLG) have attracted much more attention. Most work 
was centered on quadratic theories, i.e. on Lagrangians being quadratic polynomials 
in the Ricci tensor and the curvature scalar ||, ^, but several authors studied more 
general Lagrangians pj. These theories turned out to be inadequate as candidates for 
foundations of quantum gravity since they are non-unitary, but recently play a role as 
effective field theories. What is more relevant here, it was found that their particle 
spectrum contains a massive spin-two field. The dynamics of this field can be described 
and investigated by recasting the fourth-order NLG theory into a standard nonlinear 
second-order Lagrangian field theory. The procedure entails a decomposition of the 
dynamical data, consisting of the metric field g^i, and its derivatives up to the third 
order, into a set of independent fields describing the physical state by their values and 
their first derivatives only. In this peculiar sense, one may say that the single "unifying" 
field g^i, is replaced by (or decomposed into) a multiplet of gravitational fields. 

An adequate mathematical tool for this purpose is provided by a specific Legen- 
dre transformation |10|, [Tl| . Although the transformation has been known for more 
than a decade, it is not currently used in a systematic way. Instead, most papers on 
the nonlinear spin-two field have employed various ad hoc tricks adjusted to quadratic 
Lagrangians ||12|, |13|, |1^ (actually equivalent to the Legendre transformation for this 



particular case), but such approach does not allow one to fully exhibit the structure of 
the theory. 

Although the Legendre transformation is essentially unique, the various fields of the 
resulting multiplet can be given different physical interpretations; the different choices 
are traditionally called "frames''^. In general, the original "Jordan frame" ( JF; the name 
is borrowed from scalar-tensor theories) consisting of only the unifying metric g^^, can 
be transformed into frames including fields of definite spin in two ways. A first possibility 

^The use of the word "frame" in this sense should be deprecated, because it does not refer to the 
choice of a physical reference frame, but this abuse of terminology is now so universally adopted that 
we feel that trying to introduce here a more appropriate term, for instance "picture" as it is used in 
Quantum Mechanics, would only lead to confusion. 
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is that the field g^,y remains the spacetime metric, now carrying only two d.o.f., while the 
other degrees of freedom (previously carried by its higher derivatives) are encoded into 
auxiliary (massive) fields of definite spin: this is the Helmholtz- Jordan frame (HJF). 
Alternatively, one introduces (via an appropriate redefinition of the Legendre transfor- 
mation) a new spacetime metric g^j,,^, while the symmetric tensor is decomposed into 
spin-2 and spin-0 fields, forming in this way the massive, non-geometric components 
of the gravitational multiplet; these variables form together the "Einstein frame" (EF). 
Both frames are dynamically equivalent and very similar on the level of the field equa- 
tions: the equations of motion are second-order Lagrange field equations, and in each 
frame the corresponding spacetime metric satisfies Einstein's field equations^, thus the 
theory looks like ordinary general relativity, with the non-geometric components of the 
multiplet acting as specific matter fields. The two frames differ however in the action 
integral. In HJF the spin-0 and spin-2 fields are nonminimally coupled to gravity (to 
curvature) while there are no kinetic terms for those fields in the Lagrangian; only the 
metric has the standard Einstein-Hilbert Lagrangian R{g). In consequence, propaga- 
tion PDEs for the fields with spin zero and two arise from the action in a more involved 
way (through the variation of the metric tensor), so the theory in this frame cannot 
be obtained by minimal coupling to ordinary gravity of some additional fields already 
possessing a definite dynamics in a fixed background spacetime. Yet it is remarkable 



that in the EF variables the theory has fully standard form ||^, |T^: one recovers the 
R{g)y/—g Lagrangian for the metric and universal kinetic terms for the spin-0 and 
spin-2 fields (independently of the form of the original Lagrangian L{g) in JF), only the 
potential part of the action being affected by the actual form of L{g). The EF variables 
are uniquely characterised by these features, while in HJF different ad hoc redefinitions 
of the variables can be intertwined with the Legendre transformation (the latter being 
itself sometimes disguised as a mere change of variables) ||T2], |I3] . 

Though matematically equivalent, the two frames are physically inequivalent; the 
difference is most clearly visible while defining the energy since the latter is very sensitive 



to redefinitions of the spacetime metric |T5|. Both mathematical similarity to ordinary 



general relativity and physical arguments indicate that the EF is physical [jl5|: in HJF 
the energy-momentum tensor is unphysical, being linear in both non-geometric fields, 
while in EF the stress-energy tensor for the scalar field has the standard form and that 
for the spin-two field seems also more acceptable than in HJF. 

Anyhow, in both frames NLG theories provide a consistent description of a self- 
gravitating massive spin-two field. The field is necessarily nonlinear and in quantum 
theory it is ghostlike. The latter defect is inferred from the fact that in the linearized 
theory the Lagrangian of the field appears with the sign opposite to that for linearized 



Einstein gravity [|T^ . This fact is interpreted in classical theory as related to the occur- 



^Some authors seem instead to believe that Einstein equations can be obtained only after redefinition 
of the spacetime metric, i.e. in EF. 



3 



rence of excitations with negative energy for the field, and in consequence as a signal 
of instability of the theory. However, in this paper we show (an incomplete proof was 
given previously in |T^) that the ground state solution (vacuum) is classically stable 



at the linear level. This does not prove the stability of the vacuum state whenever 
nonlinear terms are taken into account: we stress, however, that all the above men- 
tioned features usually advocated as signals of instabihty, are equally derived within the 
linear approximation. The problem of energy is more subtle: in HJF the variational 
energy-momentum tensor is evidently unphysical, while in EF the Lagrangian is highly 
nonlinear (not even polynomial), and we will show in sect. 7 that the linear approxima- 
tion tells rather little about the energy density of the exact theory. Hence, the massive 
spin-two field generated by NLG theories is still worth investigating in the framework 
of classical Lagrangian field theory. 

Most previous works |jl2|, |13|, Q were centered on the particle spectrum of the the- 
ory and dealt only with the action integrals, while less attention was paid to the field 
equations and the structure of the theory. The field equations in HJF were given in 
13| , but these authors regarded the field equations in EF as extremely involved and 



thus intractable, they only studied the case where the spin-two field is assumed to be 
proportional to the metric tensor and thus can be described by a scalar function (its 
trace). In consequence, the dynamical consistency of the theory has always been taken 
for granted since the fourth-order equations in JF are consistent. 

In the present paper we systematically investigate the nonlinear spin-two field gen- 
erated by a NLG theory with a quadratic Lagrangian (|]) in the framework of classical 
field theory, employing a Legendre transformation. As is well known, a spin-two field 
may be mathematically represented by tensor variables of different rank and symmetry 
properties fl^, 0]. Any NLG theory generates in a natural way a representation of the 
field in terms of a symmetric, second rank tensor ip^^, = ipufj,, said this representation 
will be employed in this work. 

Although we study a concrete Lagrangian (this is motivated in sects. 2 and 3) we 
employ no tricks adjusted to it. The paper is self-contained and provides an (almost) 
full exposition of the subject. The first part describes the nonlinear spin-two and spin- 
zero fields in HJF. The particle content in this frame is well known. The main new 
results shown here are: 

• the field equations for the metric and the other two fields for any spacetime di- 
mension d > 4] 

• the fact that dimension c? = 4 is distinguished in that the scalar field is decoupled 
and can be easily removed from the theory (we do so in the rest of the paper since 
we are interested in the description of the spin-two field); 

• the fact (sect. 3) that the resulting equations of motion for the spin-two field 
do not generate further constraints besides the five ensuring the purely spin-2 
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character of the field; 

• the internal consistency of the theory in HJF, ensured by strong Noether conser- 
vation laws; 

• the linear stability of the unique ground state solution representing fiat spacetime 
and vanishing spin-two field, assessed by the fact that small perturbations form 
plane waves with constant amphtudes. 

In sect. 4 we show that the two possible ways to obtain the massless limit of the spin- 
two field described in the previous section yield the same result, so the massless limit is 
well-defined and leads to the propagation equations for gravitational perturbations in 
a Ricci-fiat spacetime. 

The main thrust of the paper is its second part (sect. 5 to 9), where we investigate 
the equations of motion in Einstein frame. Here most results are new. 

• A generic presentation of the Lagrangian theory for the nonlinear spin-two field 
in Einstein frame is given in sect. 5: in order to better exhibit the structure of 
the theory in this frame, we consider a Lagrangian being an arbitrary function of 
the Ricci tensor of the original metric in JF. We explicitly give the equations of 
motion for the massive field in the generic case, an expression (highly nonlinear) for 
the full energy-momentum tensor T^j^ of the field and four differential constraints 
imposed on the field by its dynamics. 

• This generic theory is then specialized in sect. 6 to the case of the particular 
Lagrangian of eq. (P^D, which in HJF ensured that the scalar field drops out. The 
previous, generic equations produce in this case a fifth, algebraic constraint, which 
together with those already found ensures that also in EF the massive field has 
five degrees of freedom and is purely spin-two. 

• A unique, linearly stable ground state solution is then found (without any sim- 
plifying assumptions); clearly it corresponds via the Legendre transformation to 
the ground state in HJF. The spin-two field is then redefined to make it vanish 
in this state (and in vacuum in general). It turns out that the consistency and 
hyperbolicity problems in this frame, investigated in sect. 7, are harder that in 
HJF and should be studied perturbatively; we study them in the linearized the- 
ory. The Lagrangian is computed in the lowest order (quadratic) approximation 
around the ground state solution to show the ghostlike character of the spin-two 
field. A detailed comparison is made with the theory of the linear spin-two field, 
and the fact that the energy-momentum tensor of the nonlinear theory is not ap- 
proximated by the energy-momentum tensor derived from the Wentzel Lagrangian 
is fully explained. 
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• Finally, contrary to the common belief that the full (nonlinear) system of equations 
of motion in EF are intractably involved, in sect. 8 we give a simple but nontrivial 
solution to them. 

Conclusions are in sect. 9, and Lagrange field equations and the energy-momentum 
tensor in EF for the redefined spin-two field as well as some other useful formulae are 
contained in Appendix. 

2 Equations of motion for the gravitational 
multiplet in Helmholtz— Jordan frame 

We will investigate dynamical structure and particle content of a nonlinear gravity 
(NLG) theory using Legendre transformation method |TD|, |TT|. The starting point is 
a (i-dimensional manifold M, d > 4 (later the dimensionality will be fixed to d = 4) 
endowed with a Lorentzian metric g^,^. The inverse contravariant metric tensor will be 
denoted by g'^'^ and g^"'gav = ^u'^ we introduce this nonstandard notation for further pur- 
poses. One need not view g^i, as a physical spacetime metric, actually whether g^y or its 
"canonically conjugate" momentum is the measurable quantity determining all space- 
time distances in physical world should be determined only after a careful examination 
of the physical content of the theory, rather than prescribed a priori. Formally g^^, 
plays both the role of a metric tensor on M and is a kind of unifying field which will be 
decomposed in a multiplet of fields with definite spins; pure gravity is described in terms 
of the fields with the metric being a component of the multiplet. In general dynam- 
ics for g^j^u is generated by a nonlinear Lagrangian density L\f^-g = f{g^u, Rai3fiu)\^—g 
where g = det(^^,y) and Ra/Sfiu is the Riemann tensor for ; / is any smooth (not 
necessarily analytic) scalar function. Except for Hilbert-Einstein and Euler-Poincare 
topological invariant densities the resulting variational Lagrange equations are of fourth 
order. The Legendre transformation technique allows one to deal with fully generic 
Lagrangians; from the physical standpoint, however, there is no need to investigate 
complicate or generic Lagrangians. Firstly, in the bosonic sector of low energy field 
theory limit of string effective action one gets in the lowest approximation the Hilbert- 
Einstein Lagrangian plus terms quadratic in the curvature tensor. Secondly, to obtain 
an explicit form of field equations and to deal with them effectively one needs to invert 
the appropriate Legendre transformation and in a generic case this amounts to solving 



nonlinear matrix equations. Hindawi, Ovrut and Waldram |18[ have given arguments 
that a generic NLG theory has eight degrees of freedom and the same particle spectrum 
as in the quadratic Lagrangian (|I]) below, the only known physical difference lies in 
the fact that in the generic case one expects multiple nontrival (i.e. different from flat 
spacetime) ground state solutions. This result can be also derived from the observation 
that after the Legendre transformation the kinetic terms in the resulting (Helmholtz) 
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Lagrangian are universal, and only the potential terms keep the trace of the original 
nonlinear Lagrangian. If the latter is a polynomial of order higher than two in the 
curvature tensor, the Legendre map is only locally invertible and this leads to multi- 
valued potentials, generating a ground state solution in each "branch" ; yet the form of 
the potential could produce additional dynamical contraints, affecting the number of 
degrees of freedom, only in non-generic cases. The physically relevant Lagrangians in 
field theory depend quadratically on generalized velocities and then conjugate momenta 
are linear functions of the velocities. For both conceptual and practical purposes it is 
then sufficient to envisage a quadratic Lagrangian 

L = R + aR" + bR^,{g)R^''{g). (1) 

In principle one should also include the term RapiivR"'^^'^ (in four dimensions it can be 
eliminated via Gauss-Bonnet theorem), but the presence of Weyl tensor causes troubles: 
although formally the Legendre transformation formalism works well there are problems 
with providing appropriate propagation equations for the conjugate momentum and with 
physical interpretation (particle content) of the field. We therefore suppress Weyl ten- 
sor in the Lagrangian. The Lagrangian cannot be purely quadratic: it is known from 
the case of restricted NLG theories (Lagrangian depends solely on the curvature scalar, 
L = f{R)) that the linear term R is essential and we will see that the same holds 



for Lagrangians explicitly depending on Ricci tensor The coefficients a and b have 
dimension [length]^; contrary to some claims in the literature there are no grounds to 
presume that they are of order (Planck length)^ unless the Lagrangian (|l]) arises from a 
more fundamental theory (e.g. string theory) where h is explicitly present. Otherwise 
in a pure gravity theory the only fundamental constants are c and G; then a and b need 
not be new fundamental constants, they are rather related to masses of the gravitational 
multiplet fields. Here we assume that the NLG theory with the Lagrangian (|l]) is an 
independent one, i.e. it inherits no features or relationships from a possible more fun- 
damental theory. 

As was mentioned in the Introduction, in the Legendre transformation one replaces 
the higher derivatives of the field ^^jy by additional fields. In this section we assume 
that the original field g^u keeps the role of the physical spacetime metric, and the self- 
gravitating spin-two field originates from the "conjugate momenta" to g^y. 

We recall that for a second-order Lagrangian such as (|lD one should properly choose 
the quantities to be taken as generalized velocities to define, via a Legendre map, con- 



jugate momenta [|TT|. One cannot, for instance, use the partial derivatives cjf^u^ap as 
generalized velocities since for covariant Lagrangians, e.g. (p, the Legendre map can- 
not be inverted: the Hessian, being the determinant of a 100 x 100 matrix, vanishes, 

det ( — -^^^z \ = 0. (2) 



dg^u,af3dg 
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General covariance indicates which hnear combinations of gi_iu,ai3 can be used as the 
velocities, i.e. with respect to which combinations the Lagrangian is regular (the Hessian 
does not vanish). Clearly this is Ricci tensor i?^,^. The explicit use of generally covariant 
quantities in this approach is supported by the Wald's theorem [|^, ^ that only a 
generally covariant theory may be a consistent theory of a spin-two field. Following 
jlT|, in order to decompose gfj,i, into fields with definite spins, one makes Legendre 



transformations of the Lagrangian (Q) with respect to the two irreducible components of 
Rfj^^: its trace R and the traceless part S^^. In terms of S^^, = R^i, — ^Rdfiu, S^^g^^ = 0, 
the Lagrangian reads 

L = R+{a + ^)R^ + bS^J^\ (3) 

one assumes ad + b ^ and b ^ 0. One then defines a scalar and a tensor canonical 
momentum via corresponding Legendre transformations: 

X + l = 'T-.^; (4) 

oR dSf^u 

it is convenient to identify dL/dR with x + 1 rather than with x alone. From (^): 

X = 2{a+ ^-)R and t:^" = 2bS^"' (5) 

hence fields x and vr^*^ are dimensionless and vr^*^ is traceless, Ti^^g^v = 0. The new 
triplet of field variables {g^u^X^'^^'^} defines the Helmholtz- Jordan Frame (HJF). 

Equations of motion for this frame arise as variational Lagrange equations from 



Helmholtz Lagrangian |21,P, 10, 11|. First one constructs the Hamiltonian 



H=^R + ^S,^-L (6) 
expressed in terms of g^j_y and the canonical momenta, it reads 

here = g^agvisT^""^] all indices are raised and lowered with the aid of g^'^ and g^y. 
Next one evaluates Helmholtz Lagrangian defined as 

r)T ~ 8T ~ 

Lh ^ -^R{~g, d~g, d'~g) + -^~S,u{~g. d~g, d'~g) - H{~g, x, vr), (8) 
oR dSf^u 

where the derivatives dL/dR and dL/dSp^p are set equal to the canonical momenta 
X + 1 and Ti^^ respectively, while the "velocities" R and Sfj,u explicitly depend on first 
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and second derivatives of Quiy. In classical mechanics for a first order Lagrangian L[q, q) 
one has 

LHiq,p,q,p) = pq - H{q,p) =pq-pq{q,p) - L{q, q{q,p)), (9) 

i.e. Lh is a scalar function on the tangent bundle to the cotangent bundle to the config- 
uration space; Lh does not depend on p. Similarly, in a field theory Lh is independent 
of partial derivatives of canonical momenta. In classical mechnics the action / Lndt 
gives rise, when varied with respect to p, to the equation q = dH/ dp, while varied with 
respect to q generates 

d / OLh \ OLh 



dt \ dq J dq 
the latter equation is equivalent to 



0; (10) 



p=-^ and Af^^-^ = (11) 
dq dt \ dq ) dq 

Thus Lh simultaneously generates both Hamilton and Lagrange equations of motion. 
In the case of NLG theories one is interested in replacing the fourth order Lagrange 
equations by the equivalent second order Hamilton ones. For the Lagrangian (0) Lh 
reads 

L„ = + Xfl + V - - i---... (12) 

One sees that the linear Hilbert-Einstein Lagrangian for the metric field is recovered. 
This means that Hamilton equations for g^^ are not just second order ones of any kind 
but exactly Einstein field equations^ Gfj,^ = T^u{g,x,T!')- The nonminimal coupling 
interaction terms X-R and n'^'^S^u will cause that T^,y will depend on second derivatives 
of X and n'^" and will contain Ricci tensor. Since Hilbert-Einstein Lagrangian for the 
metric field in general relativity is Lg{g) = ^R, then Lh = 2Lg + 2Lf, where Lf is 
the Lagrangian for the non-geometric components of gravity, x and ir^". The energy- 
momentum tensor is then 

^ ' 9T^u{9,X,Tr) = T^X^-gLf) (13) 



2 




~g\xR + rc^'Sa, - -rr4-^x' ' ^vr"/^7r, 



Explicitly the equations j^Lh = read 



Gf^uig) = T^^{g, X, 7r) = -xG^^u + X;iiv - 9f,u + ^^a/jvr"^^^,. 



^We use units SttG = c = 1, the signature is ( — h ++). We use all the conventions of [E2| 
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_^ a/3 7, \ TrCt _ 1„ \a 2~ _ 

2^ ■,af3yfiu -r {n-u)a 2 ^*^'" 8{ad + b) 

-^T^'^^T^apQi^v - ^7r^"7r„^ - ^^TT^v, (14) 

here f-a denotes the covariant derivative with respect to g^p and □/ = g^'^f-^u = f;fi'^- 
As remarked above T^j, more resembles the stress tensor for the conformally invariant 
scalar field p3| , ^ than that for ordinary matter. The equations of motion for x and 
Tf'^'^ are purely algebraic and clearly coincide with (||), 

These equations can be recast in the form of Einstein ones, 

G,M = 4-.^ - i(^rfT6)^^-- ^^^^ 



Comparison of eqs. (|T^) and (|I^) shows that for solutions there exists a simple linear 
expression for the stress tensor: 

T,.(^,X,vr) = ^vr,.-^^^^X^... (17) 

This relationship allows one to derive differential propagation equations for x a^nd tt'^'^. 
Before doing it we simplify the expression ( p^lf ) for T^j, with the aid of (|l^) by replacing 
R/iu by X a-iid tt'"^ and making use of the Bianchi identity for G^jy. The latter provides 
a first order constraint on x and n^'^ , 

the constraint is already solved with respect to x'^- Upon inserting (0) into the r.h.s. 
of (jl^) one gets 

m /- N 1 ~ a 1 Q ad + 2b 

Ti,y[g, X, tt) = X;/.!^ - 2 ^^t^'^ + ^ (a';'^)" " - 2{ad + b)b ^^ 

(- - i-""... + g^x' + . (19) 

Equating the trace T^^g^^^ computed from (|I^) to the trace of ( P^ ) and applying (^) 
one arrives at a quasilinear equation of motion for x, 

[4(rf - l)a + fix - ^ (^TT-^vr., + ^x^) - (^ - 2)x = 0. (20) 

10 



The field x is self-interacting and is coupled to vr^'^, the term TT^^TTa/B acts as a source for 
X- Eq- ( ^OD is a Klein-Gordon equation with a potential and an external source. The 
mass of X is 

= 4(d-l)a + d6 ^^^^ 

and may be of both signs depending on the parameters. 

To derive a propagation equation for vr^'^ one replaces derivatives of x iii (HH) by 
derivatives of tt^'^ with the help of (plsl) . Then T^j, depends on x via terms X'^fiu and 
X^dfiu', the latter is eliminated with the aid of eq. and the reappearing term is 
again removed employing (|T^). The resulting expression for T^^,, which contains terms 
XT^fiu and x^M^^' is set equal to the r.h.s. of eq. (pT]), then the terms Xdfiu cancel each 
other and finally one arrives at the following equation of motion for vr^i, , 

^ Ajad + b) ,1,1^., 

- _ 2)b ~ fe^'' -^TT^aTT^ + 

The triplet of gravitational fields is described by a coupled system of eqs. (0), (pO|) 
and (|2^ ) and the constraint (^). The equations (^) and (|2^ ) are quasilinear and 
contain interaction and self-interaction terms which cannot be removed for dimensions 
d>A. 



Since for generic dimension the dynamics of the fields tt^^ and x cannot be decoupled, 
one can obtain some information on the individual behaviour of each field by considering 
particular solutions in which only one of the two fields is excited. 

1. Let TT^jy = 0. Then eq. (p^ holds identically while the constraint (1T8|) implies 
X = const and eq. (pOD reduces to a quadratic equation, {d — 4)x^ + 2{d — 2)x = 0. 
One solves it separately for d > 4 and for d = 4. 

a) d > 4. There are two solutions, x = and x = —2{d — 2)/{d — A). From 
eq. (|16D they correspond to R^^ = and Rf,^ = - (^a-4){ad+b) Einstein 
space respectively. Therefore there exist two distinct ground state solutions: 
Minkowski space for x = and d-dimensional de Sitter space or anti-de 
Sitter for X < depending on the sign of ad + b. 

b) For d = 4 the equation has only one solution x = corresponding to R^i, = 
and there is a unique ground state solution g^^, = rj^^, (Minkowski metric) 
and vr^j. = X = 0. 
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2. Let X = 0. We consider solutions for d > 4. The scalar field equation (pO|) 
generates the algebraic equation tt^^tTq,^ = 0, besides tc'^'^-u = arising from ([T8|) . 
Eq. (^) gets reduced, upon employing the two equations, to a linear equation 

fivr^. + - 2/?„(^,)^7r°^ = 0. (23) 
This shows that the mass of the tensor field is = —1/6. 



2.1 The four— dimensional case 

Dimension four is clearly distinguished by the scalar field eq. (pOl). Setting d = 4, one 
finds two relevant properties: 

(i) both the self-interaction and the source in eq. (|20| ) vanish and x satisfies the linear 
Klein-Gordon equation 

2(3a + 6) □x-X = 0; (24) 
for 3a + 6 > the mass is real, 

ml = , ^ > 0. (25) 
^ 2(3a + 6) ^ ' 

We shall see below that this allows one to decouple completely the propagation 
equations for the fields vr^j, and x\ 

(ii) furthermore, for 3a + 6 = eq. (^) admits only one solution, x = 0, so for this 
particular choice of coefficients the scalar field disappears from the theory: we 
shall exploit this fact in the next section to concentrate our investigation only on 
the spin-two field. 



We stress that for d> 4 vanishing of the coefficient of in eq. (|20D does not imply that 
that the scalar drops out. In fact, the equation gives rise to an algebraic relationship 
between the scalar and tensor fields, 

. 4d{d-\) 2 U{d-\) 

which replaces a propagation equation for x- 

To prove the statements above, one starts from eq. (p2D, and for d = 4 one eliminates 
X from the interaction term applying eq. (p^): then is replaced by ir'^^-a/s with the 
help of (18). Then the equations of motion for vr^,^ read 

~ 2 11 
+ i(2a + 6)(3a + 6)7r°^,,^7r^, + ^^, f^^Tr^^^.^-lvr-V/s) = 0. (27) 
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The linear equation (|2^ ) for x and quasilinear eq. ( pTj ) for vr^,^ are decoupled. Equa- 
tions ( ^71) are linearly dependent since the trace (w.r.t. g^'^) vanishes identically, hence 
there are 9 algebraically independent equations (they satisfy also differential identities, 
see below). 



For the special solution x 
eq. (^) reduces to 



the differential constraint reads ti^'^ -^y = 0. Then 



□tt 



1 



al3 



Ab 



(28) 



and the mass of the field is the same as in d > 4, i.e. = —1/6. The masses of the 
massive components of the gravitational triplet, = [2(3a + h)\~^ and = — 
agree with the values found in the linear approximation to the quadratic Lagrangian for 
spin-0 and spin-2 fields by [16, 25 1 and The non-tachyon condition is then 6 < 
and 3a + 6 > 0. 

This condition shows that the term in the Lagrangian (P is essential. In fact, if 
this term is absent and the Lagrangian reads L = R + hR^j^R^^ , then (for dimensionality 
four) 



G 



Dy y 

X ^^x 



and 



□tt 



-{l + x)Trf^u 



2^a(;.i.)/37r"^ + ^g^^ (^x - ^7r"^7r„^ 



4X: 



fiu 



0. 



(29) 
(30) 

(31) 



The two fields have masses = 1/(26) and = —1/6 and one of them is necessarily 
a tachyon. It is worth noting that in the case of restricted NLG theories, i.e. L = f{R), 
the aR^ term is also essential in the Taylor expansion of the Lagrangian: for a > 
Minkowski space is stable as the ground state solution of the theory, while for a < it 
is classically unstable and for a = the solution may be unstable or stable |jl5 . 



Finally we show that all nine algebraically independent equations (^) are hyperbolic 
propagation equations for vr^^, and they contain no constraint equations. To this end 
one replaces vr^'^.j, by x''^ with the aid of (|T8|) and one arrives at the following equations: 



1 1 



'2Ra(pu)f3T<'°'^ 



2a + b 



4X: 



Aa + b 

In this form the hyperbolicity of all the equations is evident. 



'^Raifi'^^uf 

1.2 

2(3a + 6)^^^'^ - -^X^,u 



(32) 
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3 The spin— 2 field in the gravitational doublet in 
Helmholtz— Jordan frame 



We have seen in the previous section that the generic quadratic Lagrangian (P, subject 
only to the non-tachyon condition b < and 3a + 6 > 0, describes a triplet of gravita- 
tional fields, HJF= {g iJLv 1 Xi T^^^} ) where the nongeometric fields in the triplet represent 
(on the quantum-mechanical level) interacting particles with positive masses. It has 
long been known that the theory (|I|) has 8 degrees of freedom |2^, As it 

was first found in the linear approximation and then in the exact theory |[13| , 0, 



these degrees of freedom are carried by a massless spin-2 field (graviton, 2 degrees of 
freedom), a massive spin-2 field (5 d.o.f.) and a massive scalar field. 

We are interested in the dynamics and physical properties of the massive spin-2 
field. In this context, the scalar is undesirable and its presence only makes the system 
of the equations of motion more involved. One can get rid of the unwanted scalar by 
a proper choice of the coefficients in the original Lagrangian. As mentioned previously, 
for 3a + 6 = eq. (|2^) has only one trivial solution^ x = 0. We therefore restrict our 
further study to the special Lagrangian 

L = R + a{R^- m^.uR^") (33) 

and denote m? = (3a)~^ assuming a > 0. As it was noticed in |T^ this Lagrangian 
can be neatly expressed in terms of Weyl tensor, 

L = R+ ^{Lgb - Co^p^.C'^^n (34) 

where Lq-q = Rap^uR"^^^^ — AR^^R^^ + _R^, the Gauss-Bonnet term, is a total diver- 
gence in four dimensions. 

One formally repeats the operations of the previous section and replaces the expres- 
sions (pi) to (iraD by 

9 

R = Gm'^x, = — Y'^^,u, (35) 

2 

~ ~ TTl 

tt/^^^ = -?,x^' and 

Tf,u{g, X, 71") = X;f^u - - nvr^i. + '^°'(f,;u)a + — VT^^TT^^ - m'^X'^l^u - 

-~9,u hx + ^vr"^;„^ + ^vr"^vr„^) (37) 



*That the scalar degree of freedom disappears in this case was previously found in [061 n2l 
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respectively. The scalar field is still present in these equations. However the trace of 
eq. (|37D is T^^g^^^ = while the trace of eq. (|36| ) yields T^uQ^'^ = — 6m^x. Then x = 
and the scalar drops out from the theory. Although the scalar field vanishes one can- 
not, however, remove it from the Lagrangian in HJF unless one imposes the constraint 
R = already on the level of the initial Lagrangian in HF. It is more convenient to deal 
with Lagrangians containing no Lagrange multipliers and therefore the auxiliary non- 
dynamic (i.e. having no physical degrees of freedom) scalar x remains in the Helmholtz 
Lagrangian ([1^) which now reads 

2 

Lh = R + xR + t^^'S^, - 3m\2 + — TT^'^vr^,. (38) 

The system of field equations for the gravitational doublet HJF={^^i,, vr'^'"}, having 
together seven degrees of freedom, consists of Einstein field equations, 

2 

~ ~ Tfl 

Gf^uig) = Tf,^{9, tt) = — —7r^,y (39) 



and quasilinear hyperbolic propagation equations for vr 



2 

~ TTl ~ 

'f,u "^^VT^z. - 2/?«(^^,)^7r"^ + —g^uT^°''^T^aP = -"^^TT^z. - ^T^uig, Tt) = 0. (40) 

It should be stressed that, as is clearly seen from the method of deriving eqs. (pO|) 
and (p2D, the eqs. ( pOf ) and (|40|) are not simply the variational equations SLh/^q^'^ = 
and 6Lh/Sh^'^ = with the substitutions x = and R = 0. 

The field vr^*^ satisfies 5 constraints, T^^'^g^j.u = and t^^^ -y = 0. Notice that the field 
equations give rise to no further constraints. In fact, the trace of ( PP] ) and divergence of 



(|39D vanish identically if the constraints hold. A possible source of a further constraint 
is divergence of eq. (|40|). It may be shown by a direct calculation that if the equations 
91) and (HOl) hold throughout the spacetime and if the five constraints are satisfied 



everywhere, then divergence of eq. (BO) vanishes identically. This confirms the previous 



result |T^, ^ ^ 0, 1^ that this field has spin two without any admixture of lower 
spin fields. 

We now investigate the internal consistency of the theory based on eqs. (B^) and 



(^Of ). It is well known |T^, |30|, ^ that a linear spin-two field massive or massless, 
has inconsistent dynamics in the presence of gravitation since in a curved non-empty 
spacetime the field loses the degrees of freedom it had in flat spacetime and the five 
conditions which ensured there its purely spin-2 character, ^^^^1]^'^ = = ^'^"^u, are 
replaced by four differential constraints imposed on the field and Ricci tensor^. Here 

^For linear fields with spins higher than 2 it was long believed that there was no easy way 
to have physical fields on anything but Minkowski space; only recently a progress has been made for 



massive fields |32 
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we are dealing with the nonhnear spin-2 field tt^j^ and one expects that this field is 
consistent. This expectation follows from the dynamical equivalence of the Helmholtz 
Lagrangian ( pS] ) to the purely metric Lagrangian ( P^D of the fourth-order theory and 
the latter one (as well as any other NLG theory with a Lagrangian being any smooth 
scalar function of the curvature tensor) is known to be consistent. However the spin-2 
field theory in HJF should be a self-contained one and one should show its consistency 
without invoking its equivalent fourth-order version. 

We first notice a difference in the structures of the theories for the linear and nonlinear 
spin-two fields. For the linear field in Minkowski space one first derives (quite 
involved) Lagrange field equations and then either by employing the gauge invariance 
(for the massless field) or by taking the trace and divergence of the field equations for 
the massive one, one derives the five constraints ^^vV^^ = = ^^'^ ^p. The method does 
not work in a generic curved spacetime [0. For the nonlinear field vr^i, the tracelessness 



is a direct consequence of the tracelessness of S^y and of the field eq. ([T^), which now 
reads (in terms of Lf = ^{Lh — R) rather than of L^) 

XT 1 _ ™2 

- = 2^.. + = 0; (41) 



in the Helmholtz Lagrangian ( pSf ) it is not assumed that TT^,y has vanishing trace. After 
eliminating the scalar field one gets R = and the other four constraints tc^'^-u = and 
the algebraic equation for vr^,^ reduces to 

1 ~ TTl^ 

E,, = -iR,p + —n,y)=0, (42) 
while the Einstein field equations for g^j^^, are 

6 n 



-^Lh) = G,,, - f^^Cg, tt) = 0, (43) 



dgf" V2 

where T^^i, is given by eq. (^) for x = 0. Hence the constraints ensuring that vr^i, has 5 
degrees of freedom hold whenever the field equations hold. Clearly the system (^21) -(^) 
is equivalent to the system (|39|)-(^0D but the former is more convenient for studying the 
consistency of the equations. To this end one employs the coordinate invariance of the 
action integral ||30[. Under an infinitesimal coordinate transformation x'^ = x'^ + e^{x), 



\e'^\ <C 1, the metric and the spin-2 field vary as 5g'^'^ = 26^^'"^ and 

^tt'^^ = -L.Tc^"^ = -vr^^^e" + vr^'^ef^ + vr^"^;;, (44) 
here L is the Lie derivative. The action integral 

Sf= flfJ^g d'x (45) 
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is generally covariant, therefore it is invariant under the transformation 



6,Sf = j [f^,'" - Fx,^ - Fa/37r"/^;^ - 2E^^,f,7r''^ - 2F^^7r"^;,] ^ d^x = 0, (46) 

where F = 6Lf/6x = — 3m^x and 6x = — f^^'X.^- Thus the coordinate invariance 
implies a strong Noether conservation law 

f.u'" - Fx,, - E^pTi'^^,^ - 2E^^,p7T'^P - 2E,p7^^^,^ = 0. (47) 

Now assume that the field equations F = and E^^, = hold. Then also E^a-is = 
holds and the identity reduces to T^j^'^ = 0. This is a consistency condition for Einstein 
field equations (^31) . Divergence of T^,y should vanish due to the system of field equations 
and constraints without giving rise to further independent constraints. 
By a direct calculation one proves the following proposition: if the field equations x = 
and (^) and the five constraints 'n^^g^u = = ti^^ -^u hold throughout the spacetime, 
then the stress tensor given by eq. (|3^) is divergenceless, Tfj,^'^ = 0. This shows that the 
system (^2|)-(^) is consistent. 



Owing to the tracelessness of vr^'^ there is only one ground state solution for the 
system (|S^)-(|i(]|) |r3| , [T^ (i.e. the spacetime is maximally symmetric and ir'^" is co- 
variantly constant). This is Minkowski space, g^^i, = rj^^^ and vr^j, = 0. This state is 
linearly stable since small metric perturbations gfj,u = rj^^, + h^^, and excitations of vr^y 
around tt^i^ = can be expanded into plane waves h^i, = tx^^ = p^j^t, expi^ikax") with 
a constant wave vector k^, k^k" = — and a constant wave amplitude p^j, satisfying 



4 Massless spin— two field in HJF 

The massive spin-2 field n^i, has a finite range with the length scale m~^. According to 
the principle of physical continuity as the mass tends to zero, the long-range force 
mediated by 7r^,y should have a smooth limit and in this limit it should coincide with the 
strictly infinite-range theory. We therefore consider two cases: first the exactly massless 
theory resulting from the Helmholtz Lagrangian for m = and then the field equations 
of the previous section in the limit m — 0. 

After setting m = in eq. (|38D it is convenient to express the resulting Lagrangian in 
terms of Einstein tensor G^j, = S^u — \Rg,u and the trace of the spin-2 field, vr = ir^'^gfj^ty. 
One gets 

Lh = R+{x + \'rr)R + Ti^^G,,. (48) 
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Owing to the Bianchi identity this Lagrangian is invariant under the gauge transforma- 
tion 

vr/^- ^ tt''^- = vr'^- + e'^;- + e'^;/^ ^nd x ^ x' = X - ^e";a (49) 

with arbitrary vector e^. The scalar x + is gauge invariant. We notice that the scalar 
X cannot be removed already on the level of the Lagrangian since it would break 
the gauge invariance. Moreover the term xR in Lh is essential to obtain appropriate 
field equations. Hence in this approach the scalar cannot be eliminated by using the 
first principles. One can only set x = in a specific gauge. 
The field equations are 

^ = ^ = and ^ = R^,- -Ji~g^, = 0, 

which imply i?^,^ = 0. It is clear that unlike in the massive case now one cannot recover 
the original fourth-order Lagrangian (||) since the two fields are here unrelated to Ricci 
tensor. The fields x do not interact with the metric which acts as an 

empty-spacetime background metric. This suggests in turn that the two components of 
the gravitational triplet are test fields on the metric background, e.g. some excitations, 
and one may expect that they satisfy linear equations of motion. Variation of with 
respect to the metric yields Einstein field equations which are reduced in comparison to 



eq. (0), 

Gf,u = T^uig, X, 71") = X;fiu - 9fiu °X - ^ Htt^z. - ^7r"^.„^^^^ + 7r"(^.^)„ = 0; (50) 

the last equality follows from = 0. The scalar field x does not appear here in the 
combination x + ^tt as in the Helmholtz Lagrangian since upon varying with respect to 
the metric g'^'^ one has 6x = while 6tt = Ti^^5g^v {t^^^ is a fundamental field, i.e. is 
independent of the metric). 

The harmonic gauge condition tt^^ = is most convenient and the field equations 
for X and ti^^ simplify (upon employing R^^, = 0) to 

G^^y = X;f,u - 9t,u °X - ^ + Ra(,,u)pTi"^ = 0; (51) 



Their trace provides an equation for x. 



□(6x + vr) = 0. (52) 



^Another possible decomposition of g^i, in this frame into fields with definite spin, as is done in |Tj 
yields a different gauge transformation not affecting the scalar field. 
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The harmonic gauge condition does not fix the gauge uniquely and the remaining gauge 
freedom is generated by any vector satisfying the equation 

□£m + ^%m = 0- (53) 

The scalar function 6% + tt is not gauge invariant and under the transformation (^) 
varies as 6x' + 7r' = 6x + 7r — e";^- For the residual gauge freedom the scalar e".^ satisfies 
the equation Oe'^.a = which follows from (|53D upon taking its divergence. Since both 
6x + and e'^-a are solutions to the scalar wave equation, one can choose e'^ such that 
e°'.a = 6x + vr. The residual gauge freedom allows then one to use a specific gauge 
wherein 6% + tt = 0. Replacing x by tt one gets equations of motion for n^^, 

□vr^, - 2^,(^,)^7r'^^ + ^71.^, - ^g^, Dtt = 0. (54) 

The trace of these equations vanishes identically, hence there are 9 algebraically in- 
dependent equations. The scalar n is still present in (|5^ and 7r^,y is subject to four 
constraints, thus it represents 6 degrees of freedom. 

In this gauge one can take a special solution where the nondynamic field x is zero. 
Then vr = too and the field vr^i, is subject to 

□tt^, - 2i?,(^,)^7r°^ = (55) 

and the five constraints tc'^'^-i, = = tt, hence it has a definite spin equal two. 
It may be shown that eqs. (0) do not generate further constraints since their divergence 
vanishes identically provided Rf^i, = and n^'^.^, = 0. 

Secondly, one takes the limit m — in the equations of motion (^) and 
the scalar x has already been eliminated, the resulting field equations, R^u = and 
(|55D , represent the full set of solutions rather than a special class; the five constraints 
hold. 

In both cases the linear spin-2 field vr'^'^ is dynamically equivalent to gravitational 
perturbations of the background empty-spacetime metric, R^uid) = 0. If the perturbed 
metric is g^v + h^j,^, then on the level of the field equations and in the traceless-transversal 
gauge one can identify h^j_y with tt^j,. In this sense the theory for a massless spin-2 field 
TT^iy resulting from the Lagrangian (|33D is trivial. Yet fl^ claim that from the form of 
their Lagrangian for the field, which is very similar to (^S]), one can infer that the field 
is ghostlike. Since both Lagrangians in HJF have no kinetic terms, this conclusion can 
be reliably derived only in Einstein frame and this will be done in Sect. 7. 
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5 Equations of motion for the gravitational doublet 
in Einstein frame: a generic Lagrangian 

We have seen that the multiplet of fields forming HJF, in which the unifying field ^^i, 
was decomposed, is suitable for investigating the dynamical evolution, consistency and 
particle content of the theory. This, however, does not imply that the fields of this 
frame are truly physical variables. The problem of physical variables first appeared for 
a system of a scalar field interacting with gravity, since in this system one can make 
arbitrary redefinitions (see |15| and references therein). Following [|15| we say that dif- 



ferent formulations (employing different frames) of a theory provide various versions of 
the same theory if the frames are dynamically equivalent. Dynamical equivalence means 
that Lagrange equations of motion in different frames are equivalent while their action 
integral are, in general, unrelated; in classical field theory this is sufficient to regard these 
frames as various manifestations of one theory. Recently is has been shown that the two 
most important frames, JF and EF, are equivalent not only in the above classical sense: 
in the quantum context the path integral for the Lagrangian (|I|) in JF is equal to the 
path integral in EF for Einstein gravity coupled to a massive spin-2 field and a massive 
scalar |T^. In both frames many physical quantities are the same, e.g. for black holes 



all the thermodynamical dynamical variables do not alter under a suitable Legendre 
map [0 and the Zeroth Law and the Second Law of black hole thermodynamics are 
proved in EF (it is worth noticing that the proof works provided the coefficient of 
the term in the Lagrangian in JF is positive). It is not quite clear whether the latter 



proof and many other computations mentioned in |jT5[ mean that EF is the physical 
frame or merely show that it is computationally advantageous. The physical frame is 
distinguished among all possible dynamically equivalent ones in a rather subtle way: its 
field variables are operationally measurable and field excitations above the stable ground 
state solution have positive energy density (or satisfy the dominant energy condition). 
In the case of a scalar field appearing in scalar-tensor gravity theories or arising in the 
restricted NLG theories this criterion uniquely points to the physical variables (frame): 
the physical spacetime metric is conformally related (by a degenerate Legendre trans- 
formation) to the metric field of Jordan frame in which the theory has been originally 



formulated ||T^, and coincides with the EF metric. We emphasize that energy density 



is very sensitive to field redefinitions and thus is a good indicator of which variables are 
physical. In this sense the energy-momentum tensor for the spin-2 field, eq. (|39D , being 
just proportional to n^,^, is unphysicalQ. In fact, in absence of any empirical evidence 
regarding energy density for spin-2 particles, one may invoke analogy with a scalar field. 

Though a scalar field also has not been observed yet, it is generally accepted that 



^Also a different deffnitioii of a tensor tt^j/, representing the spin-two field in HJF, given in jlJ 



results in a linear energy-momentum tensor. 
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its energy-momentum tensor should be purely quadratic in the field derivatives and 
no linear terms may appear (in the kinetic part). For any long-range scalar field the 
presence of such linear terms would cause difficulties in determining its total energy. 
Given any test scalar field $ in an empty spacetime {R^y{g) = 0) one can take the 
tensor 6^^^, = — (7^,^n$, which is trivially conserved and dominates in the effective 
energy-momentum tensor at large distances and hence affects the total energy |3^. The 
degenerate Legendre transformation to the Einstein frame in scalar-tensor gravity and 
restricted NLG theories removes all linear terms and provides a fully acceptable expres- 
sion for the scalar field energy density. 



We therefore study now the other, more sophisticated way of decomposing g^i^ into a 
multiplet of fields by means of a Legendre transformation. We regard this transformation 
as a transition to the physical frame. We shall see that this transformation does not 
guarantee that the spin-2 field has positive energy density. Maybe the criterion for 
physical variables should be relaxed in the case of this field. In any case the new 
variables seem to be more physical than those in HJF. 

One generates the physical spacetime metric g^i, from gf^^ while the latter is viewed as 
a non-geometric component of gravity. It is useful and instructive to take at the outset 
the fully generic Lagrangian L = f{g^^,Ra/3) and only later to specify it to the form 
According to and [jlO| the true spacetime metric is defined asQ 



9"" = (- 



1/2 


detf 




9_ 


1/2 g^ 




KdRapJ 


dR^u 


9 


dR^u 



y yau 



(56) 



(^ff. To view 



where g = det{g^u), g = det (g^u) and g^^ is the inverse of g 
g^^ as a spacetime metric one assumes that dei{df / dRap) 7^ 0. From now on all tensor 
indices will be raised and lowered with the aid of this metric. At this point, to make the 
following equations more readable, we alter our notation and denote the original tensor 
field g^y by ip^^ and its inverse g^^ by 7^^^^. The Legendre transformation (|56D is a map 
of the metric manifold (M, t/'^j,) to another one, [M^g^^). We will not consider here the 
hard problem of whether the transformation is globally invertible^; we assume that the 
map is regular in some neighbourhood of a ground state solution. The fields g^v and 
ip^v will be referred to as Einstein frame, EF ={g^y,%l)^y}. For the generic Lagrangian 
ip^y is actually a mixture of fields carrying spin two and zero. Notice that for / as in 

t 

9_ 

®In [|l2j, and the physical metric is constructed in a more involved and tricky way; the 
outcome is equivalent to the Legendre map (p6|). 

^Some basic considerations of the problem can be found in 03]. 
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1/2 



[{1 + 2aR)-f^"' + 2bR^"'], 



(57) 
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hence for ijj^^ being Lorentzian and close to Minkowski metric, gfj_i, is also Lorentzian and 
close to flat metric (and thus invertible). This shows the importance of the linear term 
in (|I|). Furthermore, in the limit —>■ oo, i.e. 6^0, the spin-2 field vr^j, becomes a 
non-dynamic one since it is determined in terms of x and ip^i, by an algebraic equation, 

(1 + x)'^f^u + TT^aVTi." - ^^^,.vr"^7r„^ = 0. (58) 



(t^imu = in the linear approximation.) Then ( pf)) is reduced to the conformal rescaling 
of the metric (for more references cf. [|l^) being the degenerate Legendre transfor- 
mation 1^, |T^ , 

9,iu = ^lim tp^u = {I + 2aR)i)^^. (59) 

As mentioned above, for a restricted NLG theory, L = f{R), Einstein frame is the phys- 
ical one. 

For the time being we investigate the generic Lagrangian. We shall use a tensor 
being the difference of the two Christoffel connections, [T^ , 



QUa^ ^) = r°,(^) - T'^.M = {i^P,;u + - i^,u;p) , (60) 

hereafter Vq,T = Tq, denotes the covariant derivative of any T with respect to the phys- 
ical metric g^y. For any two metric tensors (not necessarily related by a transformation) 
the following identity holds for their curvatures 

This "curvature difference tensor" for Q satisfies K"^anv{4^^ fi') = and generates a "Ricci 
difference tensor" 

KM ^ R,.{i^) - R,M = Q%;a - Ql,;u + Q%Q% - QI^QL = K., (62) 
since Q^^.^y = Qau^fi^ a "curvature difference scalar" 

g^^'K,y{^P, g) = V. {g'^-^Q^^, - g'^^Ql^) + g^^ {Q%Qip - Q^^QL) ■ (63) 



Second order field equations in EF are generated by a Helmholtz Lagrangian [^, 
pO|] . First one inverts the relationship (0), i.e. solves it with respect to Ricci tensor, 
Rfj.u{4'ai3,d'ipai3,d'^4'ai3) = ^ fj,^ {g°'^ , ip a/s) ■ The fuuctious r^i, do not contain derivatives of 
g"'^ and '4'ai3- A unique solution exists providing that the Hessian 

\dRal3Rfj.u J 
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This means that / must exphcitly depend (at least quadratically) on i?^,^ and not only 
on R. It is here that the assumption that the original Lagrangian is at most quadratic 
is of practical importance. Next one constructs the Hamiltonian density (it is more 
convenient to use at this point scalar densities than pure scalars) 



and then a Helmholtz Lagrangian density, 

From (IID one finds R^.uii') = RfiAg) + K^^{Q), then 

, 1/2 



LH{g, ^) = R{g) + g^'K^^iQ) - g^^r^.ig, ^) + 



/(^/.i.,ra/3(^,^/')) 



(64) 



(65) 



(66) 



It is remarkable 



Toll that in EF the Hilbert-Einstein Lagrangian for the spacetime 
metric is recovered and the kinetic part of the Lagrangian for ^/'^jy is universal, i.e. 
is independent of the choice of the scalar function /. The only reminiscence of the 
original Lagrangian L in JF is contained in the potential part of Lh- It was far from 
being obvious that it is possible to define the physical metric in such a way that the 
gravitational part of Lh is exactly R{g). In this sense Einstein general relativity is a 
universal Hamiltonian image (under the Legendre map) of any NLG theory. The total 
divergence appearing in eq. (^) may be discarded and the kinetic Lagrangian for t/'^,^ 
reads 

K(g)^^^^(g^,gf^-g^^Qi), (67) 

what ensures that Lagrange equations of motion will be of second order. The proof that 
Ljy in EF is dynamically equivalent to L in JF is given in |T0[ and [0. As in HJF we 
write Lh = '2Lg + 2Lm with 



1 

2' 



K - ^g'^-'r,, + 



1/2 



f{lpa(3,raf3{g,i^)) 



(68) 



being the Lagrangian for ipf^^. This Lagrangian was found in a different way in |]T^ and 
14| for a different decomposition of the original metric ipfj,,^. 



Equations of motion for the metric are Einstein ones, G^y{g) 
usual 

2 6 



T^.uig.ip) where as 



-gL 



M) 



^5g^- 

However finding out T^j, directly from Lm requires a very long computation. To compute 
the dependence of T^y on the dynamical variables for solutions one can use an alternative 
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derivation. The fundamental differential relation between the two metrics, R^v{'4)a(i) = 
^fiu{g°'^ j'ipais), which represents the inverse Legendre map and is recovered in EF as one 
of the Euler-Lagrange equations, is inserted in eg. (p^) giving rise to R^v{g) = r^^ — K^y. 
Next, one takes into account that for solutions the variational energy-momentum tensor 
(or stress tensor) T^^, should be equal to the Einstein tensor, 

T^,{g, ^) = G^,{g) = -K^,{Q) + ]^g^,g''^ {K^p - r^p) + r^,. (69) 

The second equality becomes an identity upon insertion of the Legendre map. It may be 
verified that this expression coincides with that arising from the definition of the stress 
tensor, hence it is valid not only for solutions. As a consequence, the Einstein field equa- 
tions turn into an identity when g^^ and ip^^ are related by the Legendre transformation 
(|56D. From the Lagrange equations of motion for ip^j^ (see below) it will be evident that 



the use of the equations does not simplify the formula for T^^. Since Q'^^ contains 7^^^^, 
the energy-momentum tensor is a highly nonlinear function of ip^y'., in general it also 
comprises linear terms. The presence of second order derivatives of ifj^j^y in T^^, means 
that the energy density is not determined by initial data on a Cauchy surface. This is a 
generic feature of the stress tensors for vector and tensor fields: the vector gauge fields 
(represented by one-forms) and scalar fields are the only exceptions. 

The two metrics are subject to Bianchi identities. That for g^^,, G^" = 0, is of 
less practical use while the other one, V^G'^ii') = 0, may be reformulated by setting 
Rfiu = rfj,u{g,4') and replacing Vi. by Vu according to V^A" = V^A" + Q'^^A'^ for any 
vector field. In this way the identity is transformed into first order equations for r^j,. 



7"^(Va;/3 - \r^^;, " QlpT = 0. (70) 



These constitute four differential constraints on ip^u] they are independent of the equa- 
tions of motion for the field. 

The Lagrange equations of motion are derived in the standard way in a simple but 
long and strenuous computation. These read 

PM'^ _ yj ( d{2LM) \ d{2LM) _ ^_af3 ( jX^j, v),^{t^;u) ^uu^cj 

+\9'''Qip {r^Qlr + r^QZr) - P'" ^ M^'^iy^ij) - p^'' = 0, (71) 



24 



where M'^'^{'Vip) denotes the kinetic part of the equations comprising all derivative terms 
and P^'^ is the potential part, 



d 



HP 



1/2 



(72) 



The kinetic part is universal while P^'^ explicitly depends on / and rap- At first sight 
these equations seem intractably complicated but we shall see that for the special La- 
grangian (^) they may be simplified and some interesting information can be extracted 
as well as a simple solution can be found. 



Since in general ip^^y is a mixture of spin-0 and spin-2 fields, there are only four 
constraints (|TD|). Only after extracting and removing the scalar field one can derive 
a fifth constraint. In the case of the linear massive spin-2 field in flat spacetime the 
fifth constraint is generated by the trace of the equations of motion |[T^ and something 
analogous occurs in the present case though the procedure is more involved. The trace 
E^^g^v = is a second order equation while that with respect to ip^j_i, is (after many 
manipulations) 

E^'^iIj^, = r^iJ^^u - J^i^ua.^r'') - P'^'i^,. = 0, (73) 

where = ipfia'" — \'Tl"'^'^ai3;^i and r = g'^^'>p^u. We will see in the next section that 
vanishes in the particular case (^), inducing an additional constraint. 



6 Constraints and the ground state solution in EF 

From now on we restrict our study to the special quadratic Lagrangian (|33|) and expect 
that ip^y will turn out to be a massive purely spin-2 non-geometric component of the 
gravitational doublet. The Legendre transformation ( pTf ) takes the form 



with 



A = 



t 
9 

Inverting this transformation one gets 



3m^ 

1/2 



detirp)\ 



a^\l/2 



(74) 



m 
2A 



[(r - 3A)7/>^, - 



and 



R 



3m2 
~2A 



(r-4A). 



(75) 
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The Lagrangian (|63) reads now 



1 ?77^ 

Lm = -K{Q) + — (-r^ + V^'^-V'm- + 6Ar - 12A^). (76) 
Its potential part generates the potential piece of the equations of motion ([7l|) , 

7^2 _ _ i2A2\ ^ 4^^.^ + 4 (3A - r) ^^^1 (77) 

with the trace P'^'^ip^u = \m^{T - AA) = AY^Vf,^. 

Inserting r^i, from eq. (^) into the constraints ( |70D and employing an identity valid 
for any two nonsingular symmetric tensors ip^i, and g^i,, 

A,, = ^A^f^^p,^, (78) 

one arrives at = for appearing in eq. ([73|) . The constraints = 0, which are 
equivalent to (|70D, hold independently of the equations of motion. Assuming that the 
equations E^'^ = hold one gets from ([73|) that E'^'^ip^i, = = —P'^'^ip^u and this implies 
that T = AA or R{tp) = 0. Vanishing of R is evident from Einstein field equations (p9|) 
in HJF. Next inserting A = t/A in (|78D and simplifying the resulting equation by means 
of = one gets four differential constraints 



r 



2^^a° = 0. (79) 



Together with the algebraic constraint r = AA they form five constraints imposed on 
ip^u ensuring that the field has exactly five degrees of freedom and carries spin two. As 
in HJF, the gravitational doublet consists of spin-two fields, a massless (the metric) and 
a massive one |T2|, |^ |13|, |I4|, clearly its mass is the same in both frames. 



and 



The algebraic constraint reduces (|75|) and ([77| ) to 

P'" = ^ii]^' - r^MY" + W - rg^"] (80) 
It a 

2 
TTX 

Ti^u = 7^(^^A.. - 4V^A.aO- (81) 

zr 

Now one can find a simple relationship between the massive fields in both the frames. 
tt'^'^ as a function of the EF variables is (from ( [71D and (P^D) 

Tr^ng.fs.i^ap) = I det(^^)|-V2^M^ - ^M-, (82) 
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we stress that '-/^'^ is the inverse matrix to ip^i^ and all indices (including these at ip/iu) 
are shifted with the aid of g^^j. 

By definition ip^i, is a nonsingular matrix and as such it is not the most suitable 
description of a classical field which should vanish in a ground state. A field redefinition 
is required. To this aim we first determine a ground state solution. In fiat spacetime 
a ground state solution (vacuum) is Lorentz invariant, in a curvd one it should be 
covariantly constant, ipai3;fi = 0. This implies 

Q^li^ = 0' = const and i>api>'^^ = const. (83) 

The equations E^^ = reduce for this solution to P^'^ = and these read 

4^^«C = rV^M. + ii^a^r'' - \r'')g^.- (84) 

By inserting (|^) into (|8ll) one finds 



r^. = '^{r'- ^^a(^r^)g^.u = Cg^, (85) 



5r 



and the energy-momentum tensor is reduced to its potential part equal to T^j^ = —Cg^y. 
The ground state spacetime satisfies G^i, = —Cg^y and should be maximally symmetric, 
i.e. Minkowski, dS or AdS depending on the sign of C . Usually for a classical field its 
stress T^jy vanishes in the ground state; for ip^^y this amounts to = Aipafi'^"'^ and 
r^^ = 0. Multiplying (|T]) by Y'' one gets = ^rg^^ and g,,^ = 7]^,^. Then A = (r/4)2 
and the constraint r = AA yields r = 4 as r and A cannot vanish. Finally the ground 
state solution is g^^, = rj^y = ip^y. 

Next assume that C 7^ and the spacetime is dS or AdS. Any covariantly constant ip^j,^ 
can be written as ip^j^u = \Tg^iu + <PfMu with r = const, g^^(t>^u = and (t)^u;a = 0. Inserting 
the Riemann tensor for these two spacetimes into the identity R'^ a^iu4>ai3 + R'^ |3^lu4>aa = 
which follows from Ricci identity for a covariantly constant one finds = — 
dS and AdS do not admit a traceless covariantly constant Thus we have shown 
that the theory in EF has a unique ground state solution g^^i, = rj^y = ip^^. This 



result was previously found in ||T^ under the simplifying assumption that the spin-two 
field is proportional to the metric gf^j,; the proof here is generic. Of course there are 
many Einstein spaces in which = 0, e.g. setting ■j/'^y = g^i, one gets R^u{g) = 

and conversely, for 4'fj,i';a = and Rf^^ = the only solution is ipfj.u = dfiu- However 
these spacetimes are not maximally symmetric and cannot be regarded as ground state 
solutions. 

Hindawi et al. |TB] have shown stability of the ground state solution against restricted 
scalar perturbations of the field ipfiu- However to prove stability of the solution one 
should show either that arbitrary tensor excitations of the field do not grow in time. 
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or equivalently that these generic fluctuations have positive energy density. The second 
way is impractical taking into account the complexity of the stress tensor (^): we will 
return to the problem of energy density for small excitations in a forthcoming paper. The 
ground state solution is mapped by the Legendre transformation (|56D onto the ip^^, = rj^^ 
and vr^*^ = solution in HJF (sect. 3). The plane-wave perturbations of the solution in 
HJF can be transformed with the aid of the map (^21) onto similar perturbations in EF 
showing the linear stability of the vacuum in the latter frame. 

The massive spin-2 field should be described as an excitation above the ground state. 
In order to have a covariant description one makes a redefinition ip^y = g^^ + and 
assumes that the gravitational doublet is EF={g^^, In terms of one has 

furthermore r = $ + 4 and the constraints read 

$ + 4 = 4| det(5^ + <I>^)|l/^ <l>^„^" - h^, = 0. (86) 
The field equations explicitly expressed in terms of are given in Appendix. 

7 Internal structure of the theory in EF 

Consistency of the field equations in EF follows from the consistency in HJF and dy- 
namical equivalence of the two frames. Alternatively, one can prove it directly in EF 
using, as previously, the coordinate invariance. Now the proof is slightly different. Un- 
der an infinitesimal transformation x'^ = + e'^ one finds 5$^,^ = '^f]{^i-v) + '^^a^iiu'^i 
where r]^ = and = i($°^;^ + - $^.^'"). Using the definition (O) in 

6 J (M^x^—gLM = one arrives at an analogous strong Noether conservation law, 

T^"".^ - E'"^.^^^ - n^p^E""^ = 0. (87) 

The energy-momentum tensor has a purely geometric origin, i.e. setting R^^ = r^^ in 
(|69[) one gets T^^, = G^y, thus its divergence always vanishes, T^^.^y = 0. Consistency 
then requires that E^^.^p = identically if E'"^ = and the five constraints hold in 
the spacetime. However eqs. (fflD are too complicated to allow for a direct check of the 
identity. This will be done in a subsequent paper in a perturbative analysis about the 
ground state solution up to the second order. 

Also the problem of whether the Lagrange equations for are all hyperbolic prop- 
agation ones is harder in this frame. Here one should distinguish between eqs. ( [7lD and 

^"Tomboulis [Q employs a different decomposition of into and the massive spin-two field. 
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(|102| ) since the latter arise from the former upon using the constraints. After long com- 
putations one finds that all eqs. ( [7lD contain second time derivatives $^,^,00! this is due 
to the nonlinearities. It is difficult to establish whether these equations are hyperbolic. 
We will return to the problem in a subsequent paper where it will be shown that up 
to the second order in an perturbation expansion these form a nondegenerate system of 
hyperbolic propagation equations (of Klein-Gordon type). 

At present we investigate the structure of the linearized theory. Let a weak field 
be the source of some metric g^^. One should not a priori assume that g^^ is a fixed 
spacetime background for a small excitation since T^u{g, $) may contain linear terms 
and the metric will then be affected by the excitations. The linearized form of eqs. ([7l|) 
reads 



+2i?^„,^$°^ + - ($^, - ^g^,) = 0, (88) 



here we have used A ^ H-|$. The algebraic constraint gives in the linear approximation 
$ = and the differential ones reduce to ^^j/''^ = 0. Applying these constraints one 
gets the linearized version of eqs. (|102|) which one denotes by eJ^^^ = 0. It turns out 



that the linearized energy-momentum tensor (|101|) is = eJ^^^, thus for the linearized 

equations of motion (|101| - |102|) one gets = and R^u{g) = (thus is actually 

decoupled from the metric which becomes a background) and finally -E^^, is simplified 
to 

2eL = 0$^^ - + 2i?^,,^<l>"^ = 0; (89) 

clearly Ej^J coincides with upon employing the constraints. Eq. (^) is identical to 
the linearized form of eq. (^OD for n'^'^ in HJF upon replacing g^^i, by ipfj^^- This follows 
from eq. (|^) where one puts A ^ 1 and 'y^'^ ^ g^^ — . Then ^^^^ ^ ir^" and in this 
approximation one may replace by V^. 

Now one should show that the constraints $ = = ^^u'" (which for the moment 
will be referred to as secondary constraints) are preserved in time by the linear equa- 
tions (p9| ) and R^u^g) = 0. To this end one determines which equations of the system 
(|88D are not propagation ones. Using the Gauss normal coordinates in which goo = — 1 
and goi = 0,2 = 1,2, 3, one finds that the four equations Eq^J = (and equivalently 
^(1)0^4 _ g E'>^^^ = 0) do not contain the time derivatives $^,^,00! these will be re- 



ferred to as primary constraints. Hence eqs. (|88|) form a degenerate system consisting 
of 6 propagation equations E^^ = and the four primary constraints. This allows one 
to prove a proposition: 
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If the propagation equations and R^uig) = hold throughout the spacetime 



and the following constraints restrict initial data on a given Cauchy surface: 

• the four primary constraints Eq^J = 0, 

• the five secondary constraints $ = = $^1/'*^ and 

• additionally = 0, 

then all the constraints, $ = = and Eq^J = 0, are preserved in time and the 



eqs. are equivalent to 



Proof: Let the Cauchy surface have a local equation x° = in Gauss normal coor- 
dinates. The additional condition $ q = at x*^ = is essential since there are less 
primary constraints than the secondary ones. 

a) The constraint $ = 0. A propagation equation for $ is provided by taking the trace 
of eq. (P^, 2g'^^E^^ = □$ — = since Ra/3 = 0. The unique solution satisfying 
the initial conditions is $ = 0. 

b) A propagation equation for the vector = ^^^''^ arises by taking divergence of 
eq. (|89|). In fact, in empty {Rnu = 0) spacetime the contraction of the full Bianchi iden- 
tity yields Rapi^u'" = 0, then V'n^^^ = □5'^ + 2i?^^„^<l>'^^'" and 2^^ = is reduced 
to the vector Klein-Gordon equation, OS^ — m^S^ = 0. Next applying $ = one reex- 

presses the primary constraints in ( pS]) in terms of 5*^, 2i?oj^'' = Qof^Sa'"' — Sq-^— S ^-q = at 
Xq = 0. In Gauss normal coordinates these read S^^o = 0. The initial data = S^^q = 

at a;o = imply then that 5*^ = for all times. This in turn implies vanishing of Eq^^^ 
in spacetime, furthermore eqs. ( ^Qf ) and (|88|) become equivalent everywhere. 

As a byproduct it has been shown that the trace and divergence of eqs. (|89|) vanish 
identically giving rise to no further constraints besides the five secondary ones. This 
shows that the linear theory is fully consistent. It should be stressed that this theory 
is not identical with that for the linear massive test spin-2 field in empty spacetime, 
r?! , pm , though both fields are subject to the same equations of motion, ( p9D and 



Rfiu = 0, and to the same constraints. In fact, for the linear Fierz-Pauli spin-2 field 
^^y, it is known that all possible Lagrangians are equivalent (in fiat spacetime) to the 
Wentzel Lagrangian 

1 
4 



2 

m 



{^^"^^^ - , (90) 



where \1/ = rj^'^'^ ny [37, 17, pO, HOj. In contrast, the linearized theory for the field $ 



above is not self-contained as a theory for a free spin-two field in a fixed background (in 
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particular, it has not its own Lagrangian) and arises only as a limit case of the nonlinear 
theory based on the Helmholtz Lagrangian (|66| ) and ([76D. 

The difference is most easily seen while dealing with energy density. The Wentzel 
Lagrangian generates a quadratic stress tensor Tj^ for the linear field 0, while 
in the linearized theory for one finds that the lowest order terms of the expansion 
of the stress tensor are linear and vanish identically for solutions of the linearized field 
equations, while the quadratic part in the expansion of T^,y (which is not presented here) 
does not coincide with T^. 



This fact seems surprising, but one has to keep in mind that is derived from L 



w 



by replacing the fixed metric r^^^, by a generic metric g^,y and the ordinary derivatives 
with covariant ones, then taking the variational derivative with respect to g^iy. In this 
derivation one regards the variation of the spin-two field to be zero, 5g(^^v = 0. On 
the contrary, when computing the stress tensor for the exact Lagrangian Lm (0), one 
regards the field ipf^i, as fundamental, Sgipfj,i, = 0, and since ip^j^i, = g^^, + ^^j, then 
^g^^lv = —Sgfiu- The fact that the spin-two field excitations are defined to vanish when 
ip^iu is equal to the actual spacetime metric g^i, (and not to some fixed background 
metric, which would be physically objectionable and would also make the equations 
rather cumbersome) implies that the true stress tensor will contain both linear and 
quadratic terms not related to the Wentzel Lagrangian. Furthermore, even assuming 
^g^fiu = —^dfiu while taking the variation of L^y would not be enough to recover the 
quadratic terms of the true stress tensor, because the higher order terms in the expansion 
of Lm do contribute to quadratic terms in the stress tensor. The case of the linearized 
theory for clearly shows that in general the equations of motion alone are insufficient 
for determining an energy-momentum tensor for a given field 



Finally we find the expression for the Helmholtz Lagrangian (^) and (|76D in the low- 
est order approximation around the ground state solution g^^^, = rj^^, and = (where 
both metric and massive field perturbations are now taken into account). The metric 
and massive field excitations are g^u = Vf^f + ^^m^^ "^^i^ = ^ff^f + ^^^mi^- Expand- 
ing the gravitational part one finds R{g)\/—g ~ e'^Lw{h,0) (the Wentzel Lagrangian 
is in fact known to coincide with the lowest-order expansion of the Einstein-Hilbert 
Lagrangian), and similarly for 2Lm in eq. (^), so that up to a full divergence (some 
auxiliary expressions are given in Appendix) 



LniVtiu + (^hi,„,e(pf,i, + e ^i„y)^J-det{r]f,^ + ehf,^) ^ e [Lw{h,0) - Lwi^jtri)] . (91) 

Up to the second order the fields are decoupled (there are no interaction terms); also the 
field C,fj,u is absent. In this way we have rederived the well known fact that in the linear 
approximation the massive spin-two field y^^,^ arising from a nonlinear gravity theory is 



a ghost field (a "poltergeist" ) [1I6|, [1^, |IJ, 0, |l8l (a full family of unphysical ghost 



fields appears in the linearized theory according to a different approach in Pl|). This 



outcome is inescapable since Wald [^] gave a generic argument that in any generally 
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covariant theory of a number of consistently interacting spin-two fields at least one field 
is necessarily ghostlike. Thus consistency of gravitational interactions implies that the 
massive spin-two field produces states of negative norm in the state vector space of 
quantum theory. 



8 Exact solutions in Einstein frame 



The field equations ( |101H10^) are fairly involved, nevertheless one may seek for non- 



trivial solutions. It is well known (cf. e.g. |T^) that the corresponding fourth-order 
field equations for the metric g^^i, in JF admit solutions R^p{g) = which are trivial in 
the sense that they are vacuum solutions already in Einstein's general relativity. We 
therefore regard as nontrivial those solutions to eqs. (|101| - |102"D which after transforming 
back to JF yield nonvanishing Ricci tensor. 

A general method for seeking solutions is to investigate known classes of geometri- 
cally distinguished metrics depending on some arbitrary functions and check if they can 
satisfy Einstein field equations ( p.01| ) for a suitably chosen field If solutions are not 
precluded (see below) one may attempt to solve the entire system (|101| - |102| ). Clearly 



some simplifying assumptions regarding are indispensable. It turns out that the 
constraints (|86|) are very stringent and one should solve them for a given class of 
and g^i, before deahng with the field equations. 

A generic simplifying assumption is to represent the tensor in terms of few scalar 
or vector functions. The simplest ansatz, that <l>^jy is completely described by its trace, 
^iiu = j^dnu, is excluded by the constraints. In fact, the algebraic constraint (p6| 



IS 



then easily solved yielding $ = —4 or $ = 0. The first solution is excluded because it 
would correspond to vanishing of the original metric, g^i, = 0, while the second is trivial 
(in the above sense) since for ^^i, = one has Rfj.u{g) = Rfivig) = 0. 

A static spherically symmetric solution is not known. Then we consider the simpler 
case of a metric representing a plane-fronted gravitational wave with parallel rays (a 
pp wave) ^3], Such waves are characterized by a null covariantly constant wave 



vector (ray) k^, k^k^ = and = 0. First we assume that the spin-two field is of 
the form = 2k(^^W^^^, where W is some scalar function, and it is null in the sense 
that ^^yk" = 0. The latter condition holds if k^W^^ = 0, what implies that the gradient 
W^fj, is either null and proportional to fc^, or it is spacelike and orthogonal to A;^. If 
we consider the spacelike case, W'^W^^ > 0, we immediately get the trace $ = and 
then the algebraic constraint holds identically for any scalar W. The four differential 
constraints (^) reduce to the equation OW = 0. Then one finds that the expression of 
the stress tensor becomes 

T^u = -m^ (^(/.W^,.) + Iw'^'W^^k^k,^ , (92) 
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which is traceless. However, it is known that Ricci tensor for a pp wave is propor- 
tional to k^ky. This condition is met by T^^, as in (|92|) only if W'^^W^^ = 0, contrary to 
our assumption. 

The other possibility is = w k^k^ (which includes the special case of the previous 
ansatz where W^^ = wk^). Once again $ = and the algebraic constraint becomes an 
identity. The field is null by definition and the condition k^w^^ = arises now from 
the differential cosntraints. The inverse matrix is 7'^'' = g^'^ — while 

Q% = k'^k^^w,,) - ^w'^k^k, (93) 

implies 

= -Ik.k^Ow. (94) 

The energy-momentum tensor, 

Tf,u = ^k^k^ (uw - m^vo^ , (95) 

is admissible by the pp-wave metric; one then evaluates the Lagrange equations of 
motion ( |102| ): these reduce to one scalar equation 

{p -m^)w = 0, (96) 

hence T^^, = and R^^, = 0. We have arrived at a rather surprising result: although the 
spacetime is not empty since 7^ 0, the metric satisfies the vacuum field equations and 
the spin-two field necessarily behaves as a test matter and carries no energy. To elucidate 
it one takes into account that a pp-wave metric represents its own linear approximation 



around flat spacetime |l43| , the inverse matrix g^'^ equals its linearized version 

and the exact Einstein equations become linear and coincide with the linearized Einstein 
tensor. In this case linearity of the equations of motion and the assumption that the 
presence of the massive field preserves the pp-wave form of the metric imply that the 
interaction of the metric with its source is quite restricted. In fact, = w k^k^ means 
that the massive field propagates in the same direction as the pp wave or that their 
momenta are parallel. In this situation the equations of motion and T^^j are reduced to 
their linear parts. As we have seen in sect. 7, in the linearized theory T^^, = by virtue 
of E^^ = for any metric, not necessarily being a pp wave. 



To complete the solution one solves (^) for w. Writing the pp-wave metric as in 

[0 

ds'^ = —2H{u, X, y)du^ — 2du dv + dx^ + dy^^ (97) 

where the hypersurfaces u =const are null and are wave-fronts and k^ = —dfj,u, one 
finds ^ = and the Klein-Gordon eq. ( P^ ) is reduced to Yukawa equation on Euclidean 
plane. 
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Solutions to the latter one are well known, and employing polar coordinates r and 6 are 
expressed in terms of Bessel functions 

w{u^ r, 6) = \Ao{u)Jo{imr)+iBo{u)H^^\tmr)][Co{u)e + Do{u)] (99) 

oo 

+ [An{u)In{imr) + Bn{u)Kn{imr)] sin(n6') + Dn{u) cos(n6')] , 

n=l 

here and Kn are the modified Bessel functions, Jq and h'q'' are the usual Bessel and 
Hankel functions, while An-, Bn, Cn and Dn (n = 0, . . . oo) are arbitrary functions of the 
retarded time u. 

It should be emphasized that the solution (P3|)-(P^) is nontrivial in spite of i?^^, = 0. 
To show it one employs eq. (|7^); first its trace w.r. to ip^i, yields = 0, then the 
transformation is reversed to give 

2 

~ Tfi 

R"" = — —wk^kr (100) 

The ray vector is defined by A;^ = —df^u and then k'^ = g^^ku = "^^^k^, thus it is null in 
both metrics, g^^k^ky = ''j'^^k^ky = 0. This stems from the fact that ■i/'^jy is a pp-wave 
metric too of the form (|97D , with the determining function H{u,x,y) = H — |w. From 
(|97D one gets Rf^i, = k^kyUH ^ where Dif reduces to the 2-dimensional Laplacian as 
in (HI) and □ if = for i?^^ = 0. In Jordan frame one has 0^ = UH and ( |I00| ) 
follows. One checks directly that eq. (|100|) represents a solution of the fourth-order field 
equations. The ray vector is covariantly constant w.r. to ip^y too, V^/Cy = 0, since this 
property is independent of the metric functions H and H. 

Clearly in HJF Ricci tensor is also given by eq. ( [lOOl ). Then, according to (|53|), 
T^^iv = wk^ky and its energy-momentum tensor is equal, by virtue of (|36|), to Ricci 
tensor. At first sight it seems better than in EF because here a non-vanishing field 
carries a non-vanishing energy. However a stress tensor of the form ( |100| ) describes a 
flux of energy and momentum moving at light velocity, as e.g. for a plane monochromatic 
electromagnetic wave. This is not the case of a massive field (particle) of any spin. This 
stress tensor in HJF is misleading and in this sense the vanishing of T^^ in EF is closer 
to reality. However it should be emphasized that this bizarre situation is due to the 
peculiar choice = w k^k^. One should compare the stress tensors in both frames for 
more physical solutions. 



9 Conclusions 

A model for a massive spin-two field consistently interacting with Einstein's gravity 
is provided by a nonlinear gravity theory, in particular by its simplest version with 
a Lagrangian quadratic in the Ricci tensor. The consistency is achieved at the cost of 
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introducing the nonlinearity. The model can be formulated as a second order Lagrangian 
field theory in many mathematically equivalent versions ("frames"). Among them two 
are distinguished: HJF and EF. In both frames the respective spacetime metric satisfies 
Einstein's field equations with the spin-two field acting as a matter source: this is in 
fact a result of the Legendre transformation, not of the choice of a particular frame. 

Though the equations of motion in both frames are similar, the Lagrangian structure 
is different. While in EF the model has the standard form of classical field theory, in HJF 
the essential nonminimal coupling to curvature and the absence of a kinetic Lagrangian 
for the spin-two field indicate that the field should be viewed as a nongeometric member 
of the gravitational doublet (or triplet) rather that as ordinary matter. 

The transformation between the two frames is akin (in a generalized sense) to the 
known canonical transformation of Hamiltonian mechanics which interchanges the roles 
of particles' positions and momenta, {q,p) ^ {—p^Q)- In a similar way, one cannot 
assume that both the old and the new variables have the same physical interpretation 
as measurable quantities: in most cases one set of variables is in this sense unphysical. It 
is well known that for a quadratic Lagrangian (|l|) the total gravitational energy (in JF) 
is indefinite ^7\, ^ and in HJF the energy-momentum tensor for the massive spin-two 
field is unphysical as being purely linear. Thus HJF is unphysical in both senses, while 
at least in the case of restricted NLG theories (L = f{R)) Einstein frame is physical 
since the positive energy theorem holds. Thus energy plays a crucial role in determining 
the physical (measurable) variables because energy density is qualitatively sensitive to 
the change of the spacetime metric between HJF and EF. Furthermore Einstein frame 
is advantageous in that it is unique while HJF is not. 

It is not quite clear how harmful for classical nonlinear field theory is the fact that 
in the lowest order approximation the massive spin-two field is a ghost. The (exact) 
ground state solution is stable against linear perturbations. As regards energy, one 
cannot expect that the true expression for the energy-momentum tensor for the field 
(that in EF) is positive definite. That this is not the case follows from the existence 
of the solution found in sect. 8 for which T^,^ = 0. This energy-momentum tensor 
violates the commonly accepted postulate that T^^^ vanishes if and only if the matter 
field vanishes ||2^. Nevertheless one may still regard this model as a viable theory 



for a massive spin-two field. To do justice to the model one should study it in the 
second-order approximation. This will be done in a forthcoming paper. 
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Appendix 



From the constraints ( p6D one derives Q^^ = ($ + 4) ^ and g^^Q^^ = 0. Then 
Einstein field equations are 



1 



<l>,^$,,-^(7,,,$,,$'")($ + 4)-2 (101) 



+ - 



m 



2 r 



this energy-momentum tensor is much simpler than that for a linear massive spin-two 
field in Minkowski space pO]. The Lagrange field equations read 



a/3 



4 



($ + 4)"^ = 



2 r /I 



(102) 



where U^o^p = $„/3;a'^. 

In expanding the Helmholtz Lagrangian density Lnyf—g up to second order around 
the ground state solution the following expressions are useful, 



K{Q)^g^-^Lw{v,Q), 
where = "n^^'^'p^y and i = rj'^'^^^.u, 



and 



r^i^afs ^ 4 + 2e<^ + (2^ - 2^,^/i"^ + y^^^y;"^ 



(103) 
(104) 

(105) 
(106) 
(107) 
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